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We report the nonequilibrium dynamical phase transition (NDPT) appearing in a kinetic Ising 
spin system (ISS) subject to the joint application of a deterministic external field and the stochastic 
mutually correlated noises simultaneously. A time-dependent Ginzburg-Landau stochastic differ- 
ential equation, including an oscillating modulation and the correlated multiplicative and additive 
white noises, was addressed and the numerical solution to the relevant Fokker-Planck equation was 
presented on the basis of an average-period approach of driven field. The correlated white noises 
and the deterministic modulation induce a kind of dynamic symmetry-breaking order, analogous to 
the stochastic resonance in trend, in the kinetic ISS, and the reentrant transition has been observed 
between the dynamic disorder and order phases when the intensities of multiplicative and additive 
noises were changing. The dependencies of a dynamic order parameter Q upon the intensities of 
additive noise A and multiplicative noise M, the correlation A between two noises, and the amplitude 
of applied external field h were investigated quantitatively and visualized vividly. A brief discussion 
was given to outline the underlying mechanism of the NDPT in a kinetic ISS driven by an external 
force and correlated noises. 
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A. Introduction 

Ising spin system (ISS) is an ideal theoretical model 
to describe a variety of physical phenomena of a uniax- 
ial anisotropy system [1]. The nonequilibrium dynamic 
phase transition (NDPT) arises in a kinetic Ising spin 
system when a system is driven by both an external field 
and temperature simultaneously. Depending upon the 
combination of temperature and driving field, there ex- 
ist in a kinetic ISS two sorts of dynamical phases, i.e. 
symmetry-breaking order phase and symmetric disorder 
counterpart, and the two sorts of dynamical phases trans- 
form into each other with the variation of driving field 
and temperature [2 ~ 8]. The time-dependent Ginzburg- 
Landau (TDGL) equation offers, within an approxima- 
tion of mean-field, a more tractable continuum solution 
to the kinetic ISS than a discrete approach of Monte 
Carlo simulation does in evaluating the NDPT of a ki- 
netic ISS. 

Introducing a stochastic fluctuation within a kinetic 
deterministic system is a topic of interest and it arouses 
an extensive concern of multidisciplinary researchers in 
recent years. A large quantity of evidence suggests that 
the noise plays an important and constructive role in a 
nonlinear system [9-17], such as the phase transitions 
induced by a noise, the stochastic resonance (SR), the 
complexity of biology and so on. Among these, the most 
challenging subject is associated with the NDPT of a 
kinetic ISS driven by a deterministic modulation, e.g. 
a sinusoidal oscillation, and a stochastic white noise si- 
multaneously. Unfortunately, the studies on the NDPT 
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in a kinetic ISS have so far focused simply on a sin- 
gle deterministic external field [2-7], a deterministic field 
jointly coupled with either a sole Gaussian white noise or 
pulse source [8,18-20]. Few studies have ever touched on 
the NDPT of a kinetic ISS including a correlated dou- 
ble noise source. Zaikin et al. investigated the double 
SR of bistable system driven by an external determin- 
istic field with the addition of an uncorrelated double 
noise [13]. Denisov et al., however, studied the NDPT 
induced simply by a double multiplicative noise at a zero 
deterministic field [16]. Furthermore, Jia et al. reported 
the reentrant phenomenon appearing in a bistable sys- 
tem perturbed by correlated multiplicative noises at a 
zero deterministic field [14] as well as the SR due to the 
collaboration of both correlated noises and deterministic 
external field [15]. The above studies, although reveal- 
ing from the different aspects some individual features 
of the NDPT and the SR of a bistable system in differ- 
ent parameter spaces, have not yet given a clear account 
for the NDPT of a kinetic ISS, especially the feature in 
dynamic responding when the correlation of a multiplica- 
tive and additive noises, the driving field and the system 
temperature vary accordingly. 

We have provided an insight into the influence of cor- 
related additive and multiplicative noises on the growth 
rate of a tumor using the Logistic model in a previous 
paper [21]. Here we present some latest results of our 
study on the NDPT of a kinetic ISS. We visualize in a 
clear and vivid way the NDPT occurring in the multi- 
parameter space consisting of the correlation coefficient 
(A) between additive and multiplicative noises as well as 
their intensities (A, M,), and the amplitude of determin- 
istic driving field (ho). This paper is organized with five 
sections, including the current introductive section. In 
sec. 2, we give a brief description of a kinetic ISS and the 
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basic feature of the NDPT, the stochastic TDGL equa- 
tion with a correlated multiplicative and additive noise 
source as well as the algorithm for solving numerically the 
relevant Fokker-Planck equation. The computational re- 
sults and the relevant analysis, especially the mechanism 
underlying the peculiar reentrant NDPT in a kinetic ISS, 
are presented in detail in sec3 and 4, respectively. The 
characteristics of the NDPT in a kinetic ISS are summa- 
rized in sec. 5. 



(t(t)p(t')) = 2\VMAS(t - t') (3-3) 

Here M, A and A in eqs.(3) symbolize the intensities 
of mutually correlated multiplicative and additive noises, 
and the correlation intensity between them, respectively. 

To quantify the process of dynamic phase transition, 
the dynamic order parameter Q is defined as 



B. Nonequilibrium dynamical phase transition and 
theoretical model 



Considering the kinetic Ising model of N interacting 
spins driven by an external deterministic field, one can 
express its Hamiltonian as 



H = --J2S i -S j -h(t)J2Si 



(1) 



{hi} 



where J is the exchange-coupling constant and spin S 
= ±1. Symbol i and {i,j} stand for a single spin-site and 
a spin-pair of the nearest neighbor, respectively. The 
time variation of driving field h(t) is written in a simple 
cosinoidal form. The stochastic TDGL equation of Eq. 
(1) within the approximation of mean-field is described 
by 



Q = — d> m{t)dt 

Z7T , 



(4) 



Actually, dynamic order-parameter Q is the period- 
average of the instant order-parameter m over evolu- 
tion time t. Q=0 and Qt^O correspond to a symmetric 
dynamic disorder phase and a symmetry-breaking dy- 
namic order phase, respectively. And they were also 
referred to as the symmetry-restoring oscillation (SRO) 
and symmetry-breaking oscillation (SBO) [8]. The pur- 
pose of the current paper is to investigate quantitatively 
the dependence of Q parameter of a kinetic ISS upon M, 
A, A and ho combined. 

The relevant Fokker-Planck (F-P) equation of eq.(2) is 
given as with the initial and boundary conditions speci- 
fied 



dptrn^ = _ 2 d[C(m,t)p{ m ,t)] + P 
dt dm dm 2 - 



(5-1) 



dm ^SH(m) . . . . 

— = -r-±J. + m &) + m (2-i) 



C{m,t) = (1 + M)m-m 3 + h cos(ujt), (5-2(1)) 
B(m) = Mm 2 + 2\V~MAm + A (5-2(2)) 



H(m) is the usual m 4 Hamiltonian for the instant 
order-parameter m(t) and consists of the following terms: 



uf \ I jd r r 2 , U 

Him) — I d r\ — m H ? 

J 2 4 



m ■ h cos(wt)} (2-2) 



r and uo are two phcnomenological coefficients. The 
reduced temperature ro is defined as ^T-Tco;T and Tco 
are the temperature of system and static critical tem- 
perature, respectively; ho and to denote the amplitude 
and frequency of a deterministic driving field; t is the 
evolution time of the system; £(t) and p(t) represent the 
multiplicative and additive Gaussian noises which satisfy 
the zero-mean and autocorrelation, expressed as 

(£(*)> - 0, mm) = 2M5(t - t') (3-1) 



(p(t))=0,(p(t)p(t'))=2A5(t-t') (3-2) 



p{m,0) = S(m - 0),p(-oo,t) =p(oo,i) = (5-3) 

Apparently, partial differential eq. (5-1), owing to its 
explicit time-dependence, is non-autonomous in nature 
and it is feasible to work out the numerical solution of 
F-P eq. (5-1), i.e. probability distribution of m against 
time t, p(m,t). 

According to the definition of Q parameter in eq. (4) 
as well as the numerical solution of p(m,t) above, we can 
figure out finally the dynamic order parameter Q of a 
kinetic ISS using 



Q 



m(t)p(m,t)dtdn 



F-P 



/;/: 



(6) 



p(m,t)dtdm 



Fig. la shows an asymptotically steady and symmet- 
ric p(m,t) at a zero driving field and the zero corre- 
lation of two sorts noises. Time-dependent oscillations 
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of p(m,t) are displayed in Fig. lb (nonzero driving field 
but zero correlation of noises ) and lc (nonzero driving 
field and nonzero correlation). The symmetry of p(m,t) 
is destroyed by the mutually correlated noises and no 
asymptotically steady solution is available when a time- 
dependent driving field is involved. 



C. Results and analysis 

Fig. 2 displays the dependence of the reduced Q pa- 
rameter upon both M and A under the condition of a 
different correlation intensity A. The dynamic order pa- 
rameter Q, due to the cooperation of the driving field 
with the correlated noises, attains a distinctive resonant 
peak within a certain range of M and A, quite similar to 
the characteristic of a common stochastic resonance when 
the ratio of signal versus noise is evaluated. A larger cor- 
relation intensity leads to a higher Q peak as shown in 
Fig. 2. The existence of a peaking Q parameter indi- 
cates the enhancement of dynamic ordering in a kinetic 
ISS if the deterministic driving force could synchronize 
with the stochastic force. Judging from the viewpoint 
of the NDPT, the resonance-like dependence of Q pa- 
rameter is a kind of reentrant phenomenon that occur 
usually in a peculiar noise circumstance. Both additive 
and multiplicative noises give rise to a reentrant trend of 
the dynamic order parameter but, according to Fig. 2, 
the multiplicative noise has a more evident effect than 
the additive one does. 



Fig. 3 exhibits two four-dimensional panoramic plots 
of the sliced contour-lines of dynamic order Q within 
the parameter space of A, M, A and A, M, h ,. It has 
been displayed vividly that a correlated noise source in- 
duces the enhancement of dynamic ordering as well as the 
NDPT of reentrant characteristic. In contrast to noise 
intensity, the amplitude of a deterministic driving field 
simply results in a monotonic trend of Q parameter. 



D. Discussion 

When an interacting collaboratively many-body sys- 
tem, such as the ISS in the present paper, is driven in 
time by an external field, the system cannot respond in- 
stantaneously to the perturbation due to the relaxational 
delay of the system itself [2], and therefore there come 
along a series of intriguing noncquilibrium phenomena, 
such as dynamic response, dynamic hysteresis and dy- 
namic phase transition. There exist two competing time 
scales in a kinetic spin system, namely the oscillation 
period of an external driving field t w = 27r/wand the 



intrinsic relaxation time of a spin system t s . The lat- 
ter tends to diverge because of the critical slowing-down 
when temperature approaches to a critical point. The 
NDPT takes place only when the two competing time 
scales, t u and r s , match each other someway in quantity 
[3,4]. The involvement of a noise source introduces the 
third competing time scale tk- As for an activated Brow- 
nian particle jumping between the double potential wells, 
the noise-induced hopping between the local equilibrium 
states with the Kramers rate "f K as given as [9] , 



Ik =7o ex P(— p-) (7) 

where 7 is a prefactor, AV and D are the height of a 
potential barrier and a noise intensity, respectively. Obvi- 
ously, tk — ^/1k an d therefore tk is adjustable through 
the noise intensity D. The match between t w and tk 
brings about the common SR [9] , whereas the match be- 
tween t s and tk causes a noise-induced nonequilibrium 
phase transition. It is the match of three competing time 
scales, t s , r w and tk, that gives rise to the nonequilib- 
rium dynamic phase transition which we address in this 
paper. For the sake of simplicity, we assume a fixed am- 
plitude and frequency of the driving field and a constant 
reduced temperature. Hence, r w and t s are invariable in 
discussion below. Set u> = 2 then r w = n. According to 
either the theory of para- and ferro-magnetic resonance 
[22] or the linear responding theory [23], one can work 
out an analogic temperature dependence of t s out of ei- 
ther theory above. The t s we figure out by means of the 
linear responding theory [23] is between 2.5~6 within a 
wide range of temperature. We plot in Fig. 4 the vari- 
ation of tk versus the additive noise intensity A using 
a reduced equation 7, tk = 1/(7at/7o)- For the TDGL 
model given by eq. (2-2), AV = 1/4. It is easily observed 
from Fig. 4 that if A < 0.25, the tk is nearly divergent, 
and when 0.25 < A < 1, tk is between 1.5~10 which 
matches with r w and r s in an order of magnitude. When 
A > 1, tk = 1- The trend of tk against the additive 
noise intensity A in Fig. 4 complies excellently with the 
Q~A trend demonstrated in Fig. 2. The optimal additive 
noise Ap at which the dynamic order parameter Q has 
a peak value is around 0.25~0.5. The effect of a multi- 
plicative noise, similar to the trend of an additive noise 
as mentioned above, is more noticeable in that the hop- 
ping of a Brownian particle agitated by a multiplicative 
noise is even modulated by the instant order parameter 
m(t) itself. Because of \m(t)\ ^ 1 as well as the m(t) 
declines with the rise of temperature, the multiplicative 
noise must satisfy the requirement of Mp > Ap so as 
to achieve the same effect that an additive noise has. 
Fig. 3 shows that the optimal multiplicative noise Mp is 
about 3, greater than its counterpart Ap. Considering 
that the modulation the order parameter m impacts on 
the multiplicative noise intensity M, one can understand 
why the dynamic order parameter Q relies upon the M 
more drastically than the A when M > Mp. And that is 
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the reason why the range of reentrance induced by a mul- 
tiplicative noise is much narrower than that caused by an 
additive noise, as typically shown in Figs. 2 and 3. The 
noise-induced dynamic ordering and the NDPT come into 
being only when tk complies someway with r^and t s 
in quantity. In a word, the matches among the three 
competing time scales, what we deal with and present 
in this paper, are the fundamental condition for the oc- 
currence of a dynamic ordering, and hence the nonequi- 
librium dynamic phase transition in a kinetic Ising spin 
system subject to a deterministic external field and a 
correlated noise source simultaneously. 



additive and multiplicative noises, the relevant correla- 
tion intensity between the two sorts of noises together, 
the dynamic order parameter take an apparent reentrant 
trend against noise intensity, indicating the development 
and the fading away of dynamic ordering due to the vari- 
ation from a cooperative match to a mismatch among the 
deterministic driving field, stochastic source and the in- 
trinsic relaxation of a spin system itself. A multiplicative 
noise has a more obvious effect on inducing the dynamic 
order in a kinetic Ising spin system than an additive noise 
does. 



E. Conclusions 



The kinetic Ising spin system, when driven by a de- 
terministic external field and a stochastic source of cor- 
related additive and multiplicative noise simultaneously, 
holds a characteristic of stochastic resonance and gives 
rise to the noncquilibrium dynamic phase transition. 
Within a multi-parameter space consisting of the am- 
plitude of a deterministic external field, the intensities of 
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FIG. 1: figurel: (a) A=M=0.5, hO=0.0, A=0.0 



FIG. 2: figurel: (b) A=M=0.5, h0=0.2, A=0.0 



FIG. 3: figurel: (c) A=M=0.5, h0=0.2, A=0.5 
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FIG. 4: Fig. l.The distribution function of probability p(m,t) 
versus instant order parameter m and evolution time t at zero 
driving field and zero correlation of noise (a), nonzero driving 
field and zero correlation of noise (b), and nonzero driving 
field and nonzero correlation of noise (c). 



FIG. 5: figure2: (a) 
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FIG. 6: figurc2: (b) 



FIG. 7: figure 2: (c) 



FIG. 8: figure 2: (d) Fig2. The dependence of 

dynamic order Q upon additive noise intensity A and multi- 
plicative noise M at different correlation intensities of noises 
(A=0.2, A=0.5). 



FIG. 9: figure 3: (a) 
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FIG. 10: figure 3: (b) Fig. 3. The sliced contour 

of dynamic order parameter Q within the parameter space of 
A, M and A (above), of A, M, and hO (below) 



FIG. 11: Figure 4. Reduced tk versus additive noise intensity 
A. 
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